A relativistic quark model with infrared confinement 

and the tetraquark state 

S. Dubnicka*, A.Z. Dubnickova^ M.A. Ivanov**, J.G. Korner* and 

G.G. Saidullaeva § 

* Institute of Physics Slovak Academy of Sciences Dubravska cesta 9 SK-842 28 Bratislava, Slovak Republic 
^ Comenius University Dept. of Theoretical Physics Mlynska Dolina SK-84848 Bratislava, Slovak Republic 
** Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141 980 Dubna, Russia 
^Institut fur Physik, Johannes Gutenberg-Universitat, D-5 5099 Mainz, Germany 
^Al-Farabi Kazak National University, 480012 Almaty, Kazakhstan 

Abstract. We explore the consequences of treating the X(3872) meson as a tetraquark bound state. As dynamical framework 
we employ a relativistic constituent quark model which includes infrared confinement in an effective way. We calculate the 
decay widths of the observed channels X — > J/\[f + 2tz(3k) and X —¥ D° + D° + 71° via the intermediate off-shell states 
X — > J /\jf + p(ffl) and X — > D + D* . For reasonable values of the size parameter A^ of the X(3872) we find consistency with 
the available experimental data. 
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INTRODUCTION 

A narrow charmonium-like state X(3872) was ob- 
served in 2003 in the exclusive decay process B ± — > 
K±n + n~J/y [1]. The X(3872) decays into n+n-J/y 
and has a mass of mx = 3872.0 ±0.6(stat) ±0.5(syst) 
very close to the M D o +M D ,o = 3871.81 ±0.25 mass 
threshold [2]. Its width was found to be less than 2.3 
MeV at 90% confidence level. The state was confirmed 
in B-decays by the BaBar experiment [3] and in p~p pro- 
duction by the Tevatron experiments [4] . 

From the observation of the decay X(3872) — > J/y/y 
reported by [5], it was shown that the only quantum 
numbers compatible with the data are J PC = 1 ++ or 2 h . 
However, the observation of the decays into Z) S°^ by 
the Belle and BaBar collaborations [6] allows one to 
exclude the choice 2 h because the near-threshold decay 
X — > D°D°n is expected to be strongly suppressed for 
J = 2. 

The Belle collaboration has reported evidence for the 
decay mode X —> n + 7i^7i J /y/ with a strong three-pion 
peak between 750 MeV and the kinematic limit of 775 
MeV [5], suggesting that the process is dominated by the 
sub-threshold decay X — > coJ /y/. It was found that the 
branching ratio of this mode is almost the same as that of 
the mode X — > n + n~J /y. 
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These observations imply strong isospin violation be- 
cause the three-pion decay proceeds via an intermedi- 



ate o-meson with isospin whereas the two-pion decay 
proceeds via the intermediate p -meson with isospin 1. 
Also the two-pion decay via the intermediate p -meson is 
very difficult to explain by using an interpretation of the 
X(3872) as a simple cc charmonium state with isospin 0. 

There are several different interpretations of the 
X(3872) in the literature: a molecule bound state 

(D°D* ) with small binding energy a tetraquark state 
composed of a diquark and antidiquark, threshold cusps, 
hybrids and glueballs. A description of the current 
theoretical and experimental situation for the new 
charmonium states may be found in the reviews [7] . 

We provided in Ref. [8] an independent analysis of the 
the properties of the X(3872) meson which we interpret 
as a tetraquark state as in [9]. We worked in the frame- 
work of the relativistic constituent quark model which 
has recently been extended to include infrared confine- 
ment effects [10]. 



THEORETICAL FRAMEWORK 

The authors of [9] suggested to consider the X(3872) 
meson as a J c = 1 ++ tetraquark state with a sym- 
metric spin distribution: [cgj^o [c^]j = i + [c^]s = i [cq]s=o< 
(q = u,d). The nonlocal version of the four-quark inter- 



polating current reads 



diagram will be represented as 
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where w\ = w 2 = m c /2{m a + m c ) and wj, — W4 = 
m q /2{m q + m c ). The matrix C = y / 2 is the charge con- 
jugation matrix. The effective interaction Lagrangian de- 
scribing the coupling of the meson X q to its constituent 
quarks is written in the form 

Jfi n t = 8xXqn(x)-Jx q {x), (q = u,d). (3) 

The state X u breaks isospin symmetry maximally so the 
authors of [9] take the physical states to be a linear 
superposition of the X u and Xj states according to 

Xi = Xj OVJ = X u cos 9 + Xj sin 9 , 
X h = X high = -X u sin 9 + X d cos 9 . (4) 

The mixing angle 9 can be determined from fitting the 
ratio of branching ratios Eq. (1). 

The coupling constant gx in Eq. (3) will be deter- 
mined from the compositeness condition Zh = 0, see, 
e.g. Refs. [11] and [12]. It gives 

Z x = l-U' x (m 2 x )=0, (5) 

where Fix (p 2 ) is the scalar part of the vector- meson mass 
operator. The corresponding three-loop diagram describ- 
ing the X-meson mass operator is shown in Fig. 1. We 



c 




FIGURE 1. Diagram describing the Xu-meson mass opera- 
tor. 

will choose a simple Gaussian form for the vertex func- 
tion with the only dimensional parameter Ax character- 
izing the size of the X-meson. 

In [10] we described how to integrate n -point one-loop 
diagrams and how to implement infrared confinement of 
quarks in this process. We extend our loop integration 
techniques to the case of arbitrary number of loops. Let 
n, £ and m be the number of the propagators, loops 
and vertices, respectively. In Minkowski space the ^-loop 



H(pi,...,p m ) = 

j* m n 
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K +n = e$k, w;i) 2 (6) 
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where the vectors fc, are linear combinations of the loop 
momenta The v, are linear combinations of the ex- 
ternal momenta /?, to be specified in the following. The 
strings of Dirac matrices appearing in the calculation 
need not concern us since they do not depend on the mo- 
menta. The external momenta pi are all chosen to be in- 

m 

going such that one has £ pi = 0. 

i=l 

Using the Schwinger representation of the local quark 
propagator one has 

S(k) = (m+ k)Jdj3 e-P {m2 - k2) . (7) 


The integrand in Eq. (6) has a Gaussian form and may 
be integrated out explicitly. After doing the loop integra- 
tions one obtains 

n = y>j3FQ3i,...,j3 n ), (8) 


where F stands for the whole structure of a given dia- 
gram. The set of Schwinger parameters j3; can be turned 
into a simplex by introducing an additional f-integration 
leading to 

1 

n = jdtt"- l Jd"ad(l-Y d a^F(ta l ,...,ta„). (9) 

1=1 

There are altogether n numerical integrations: (« — 1) a- 
parameter integrations and the integration over the scale 
parameter t . The very large t -region corresponds to the 
region where the singularities of the diagram with its lo- 
cal quark propagators start appearing. However, as de- 
scribed in [10], if one introduces an infrared cut-off on 
the upper limit of the t-integration, all singularities van- 
ish because the integral is now convergent for any value 
of the set of kinematic variables. By introducing the in- 
frared cut-off one has removed all potential thresholds in 
the quark loop diagram, i.e. the quarks are never on-shell 
and are thus effectively confined. We take the cut-off pa- 
rameter A to be the same in all physical processes. 

Next we evaluate the matrix elements of the transitions 
X J/y + p (co) and X ->• D + D*. The relevant Feyn- 
man diagrams are shown in Fig. 2. Since the X(3872) is 




FIGURE 2. Feynman diagrams describing the decays X 
J/y + p{m) andX ^D + D*. 



The theoretical value for this rate ratio depends only 
weakly on the size parameter Ax 



T(X D°D° 7t ^ 



theor 



6.0±0.2. 



(12) 



The theoretical error reflects the Ax dependence of the 
ratio. The ratio lies within the experimental uncertainties 
given by Eq. (11). 



very close to the respective thresholds in both cases, the 
intermediate p, co and D* mesons have to be treated as 
off-shell particles. 



NUMERICAL ANALYSIS 

Using the calculated matrix elements for the decay X —> 
J/\j/ + p((o) one can evaluate the decay widths X — > 
J/y/+2n(37t). We employ the narrow width approxima- 
tion for this purpose. 

The adjustable parameters of our model are the con- 
stituent quark masses m q , the scale parameter A charac- 
terizing the infrared confinement and the size parameters 
Am- They were determined by using a least square fit to 
a number of physical observables, see [10]. 

There are two new free parameters: the mixing angle 
9 in Eq. (4) and the size parameter Ax- We have varied 
the parameter Ax in a large interval and found that the 
ratio 

r(X u ^J/\jf+3%) 



0.25 



(10) 



r(x u ^j/Y+27c) 

is very stable under variations of Ax- Hence, by using 
this result and the central value of the experimental data 
given in Eq. (1), one finds 9 « ±18.4° for X t ("+") 
and Xh ("-"), respectively. This is in agreement with the 
results obtained in both [9]: 9 w ±20° and [13]: 9 w 
±23.5°. The decay width is quite sensitive to the change 
of the size parameter Ax- A natural choice is to take a 
value close to Aj/y and Ajj c which are both around 3 
GeV. We have varied the size parameter Ax from 3 up to 
4 GeV and found that the decay width T{X —> J/ y+n 7t) 
decreases from 0.30 up to 0.07 MeV, monotonously. 
This result is in accordance with the experimental bound 
T(X(3872)) < 2.3 MeV and the result obtained in [9]: 
1.6 MeV. 

In a similar way we calculate the width of the decay 
X -> D°D° n° which was observed by the Belle Coll. 
and reported in [6]. We have varied Ax from 3 up to 4 
GeV and found that the decay width r(X, ->■ D°D° 7t°) 
decreases from 1.88 up to 0.41 MeV, monotonously. 

Using the results of [2], one calculates the experimen- 
tal rate ratio 
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